Search for: All records

ABSTRACT We study the morphology of hundreds of simulated central galaxies in the stellar mass range $$M_\star =$$ 107.5–1011  $$\rm M_\odot$$ from the FIREbox cosmological volume. We demonstrate that FIREbox is able to predict a wide variety of morphologies, spanning from disc-dominated objects to spheroidal galaxies supported by stellar velocity dispersion. However, the simulations predict a strong relation between morphology (degree of rotational support) and stellar mass: galaxies comparable to the Milky Way are often disc-dominated while the presence of stellar discs mostly vanishes for dwarfs with $$M_\star < 10^9 ~$$\rm M_\odot$$. This defines a ‘morphology transition’ regime for galaxies with $$10^9 < M_\star /\rm {M_\odot }< 10^{10}$$ in which discs become increasingly common, but below which discs are rare. We show that burstiness in the star formation history and the deepening of the gravitational potential strongly correlate in our simulations with this transition regime, with discs forming in objects with lower levels of burstiness in the last $$\sim 6$$ Gyr and haloes with mass $$\sim 10^{11} ~ \rm {{\rm M}_{\odot }}$$ and above. While observations support a transition towards thicker discs in the regime of dwarfs, our results are in partial disagreement with observations of at least some largely rotationally supported gas discs in dwarfs with $$M_\star < 10^9$$\rm M_\odot$$. This study highlights dwarf morphology as a fundamental benchmark for testing future galaxy formation models. 

Abstract We utilize the cosmological volume simulation FIREbox to investigate how a galaxy’s environment influences its size and dark matter content. Our study focuses on approximately 1200 galaxies (886 central and 332 satellite halos) in the low-mass regime, with stellar masses between 10⁶and 10⁹M_⊙. We analyze the size–mass relation (r₅₀–M_⋆), the inner dark matter mass–stellar mass ( $M_{\mathrm{DM}}^{50}$–M_⋆) relation, and the halo mass–stellar mass (M_halo–M_⋆) relation. At fixed stellar mass, we find that galaxies experiencing stronger tidal influences, indicated by higher Perturbation Indices (PI > 1) are generally larger and have lower halo masses relative to their counterparts with lower Perturbation Indices (PI < 1). Applying a Random Forest regression model, we show that both the environment (PI) and halo mass (M_halo) are significant predictors of a galaxy’s relative size and dark matter content. Notably, becauseM_halois also strongly affected by the environment, our findings indicate that environmental conditions not only influence galactic sizes and relative inner dark matter content directly, but also indirectly, through their impact on halo mass. Our results highlight a critical interplay between environmental factors and halo mass in shaping galaxy properties, affirming the environment as a fundamental driver in galaxy formation and evolution. 

Abstract: Atacama Large Millimeter/submillimeter Array observations have shown that candidate “post-starburst” galaxies (PSBs) at z~0.6 can retain significant molecular gas reservoirs. These results would imply that—unlike many model predictions—galaxies can shut down their star formation before their cold gas reservoirs are depleted. However, these studies inferred star formation rates (SFRs) either from [OII] line fluxes or from spectral energy distribution (SED) modeling and could have missed large dust-obscured contributions to the SFRs. In this study, we present Keck/NIRES observations of 13 massive (M_* >= 10^11M_⊙) PSBs, which allow us to estimate Hα SFRs in these gas-rich PSBs. We confirm the previously inferred low SFRs for the majority of the sample: 11/13 targets show clear Hα absorption, with minimal infilling indicating dust-corrected SFRs of <4.1Msun/yr. These SFRs are notably low given the large H2 reservoirs (∼(1–5) × 10^10Msun) present in 5/13 of these galaxies, placing them significantly offset from star-forming galaxies on the Kennicutt–Schmidt relation for star-forming galaxies. The [NII]/Hα ratios of all 13 PSBs imply contributions from non-star-forming ionization mechanisms (e.g., active galactic nuclei, shocks, or hot evolved stars) to their Hα emission, suggesting that even these low ongoing SFRs may be overestimated. These low Hα SFRs, dust corrected using Av estimates from SED fitting, confirm that these galaxies are very likely quiescent and, thus, that galaxies can quench before their cold gas reservoirs are fully depleted. 

ABSTRACT Galaxy formation is a complex problem that connects large-scale cosmology with small-scale astrophysics over cosmic time-scales. Hydrodynamical simulations are the most principled approach to model galaxy formation, but have large computational costs. Recently, emulation techniques based on convolutional neural networks (CNNs) have been proposed to predict baryonic properties directly from dark matter simulations. The advantage of these emulators is their ability to capture relevant correlations, but at a fraction of the computational cost compared to simulations. However, training basic CNNs over large redshift ranges is challenging, due to the increasing non-linear interplay between dark matter and baryons paired with the memory inefficiency of CNNs. This work introduces EMBER-2, an improved version of the EMBER (EMulating Baryonic EnRichment) framework, to simultaneously emulate multiple baryon channels including gas density, velocity, temperature, and H i density over a large redshift range, from $z=6$ to $z=0$. EMBER-2 incorporates a context-based styling network paired with Modulated Convolutions for fast, accurate, and memory efficient emulation capable of interpolating the entire redshift range with a single CNN. Although EMBER-2 uses fewer than 1/6 the number of trainable parameters than the previous version, the model improves in every tested summary metric including gas mass conservation and cross-correlation coefficients. The EMBER-2 framework builds the foundation to produce mock catalogues of field level data and derived summary statistics that can directly be incorporated in future analysis pipelines. We release the source code at the official website https://maurbe.github.io/ember2/. 

ABSTRACT Atomic hydrogen (H i) serves a crucial role in connecting galactic-scale properties such as star formation with the large-scale structure of the Universe. While recent numerical simulations have successfully matched the observed covering fraction of H i near Lyman Break Galaxies (LBGs) and in the foreground of luminous quasars at redshifts $$z \lesssim 3$$, the low-mass end remains as-of-yet unexplored in observational and computational surveys. We employ a cosmological, hydrodynamical simulation (FIREbox) supplemented with zoom-in simulations (MassiveFIRE) from the Feedback In Realistic Environments (FIRE) project to investigate the H i covering fraction of Lyman Limit Systems ($$N_{{\text{H}}\, \rm{{\small I}}} \gtrsim 10^{17.2}$$ cm$$^{-2}$$) across a wide range of redshifts ($z=0-6$) and halo masses ($$10^8-10^{13} \, \,\mathrm{ M}_{\odot }$$ at $z=0$, $$10^8-10^{11}\, \,\mathrm{ M}_{\odot }$$ at $z=6$) in the absence of feedback from active galactic nuclei. We find that the covering fraction inside haloes exhibits a strong increase with redshift, with only a weak dependence on halo mass for higher mass haloes. For massive haloes ($$M_{\mathrm{vir}} \sim 10^{11}-10^{12} \,\mathrm{ M}_{\odot }$$), the radial profiles showcase scale-invariance and remain independent of mass. The radial dependence is well captured by a fitting function. The covering fractions in our simulations are in good agreement with measurements of the covering fraction in LBGs. Our comprehensive analysis unveils a complex dependence with redshift and halo mass for haloes with $$M_{\mathrm{vir}} \lesssim 10^{10} \,\mathrm{ M}_{\odot }$$ that future observations aim to constrain, providing key insights into the physics of structure formation and gas assembly. 

ABSTRACT Our ability to trace the star-forming molecular gas is important to our understanding of the Universe. We can trace this gas using CO emission, converting the observed CO intensity into the H$$_2$$ gas mass of the region using the CO-to-H$$_2$$ conversion factor ($$X_{\rm{{\small CO}}}$$). In this paper, we use simulations to study the conversion factor and the molecular gas within galaxies. We analysed a suite of simulations of isolated disc galaxies, ranging from dwarfs to Milky Way-mass galaxies, that were run using the fire-2 subgrid models coupled to the chimes non-equilibrium chemistry solver. We use the non-equilibrium abundances from the simulations, and we also compare to results using abundances assuming equilibrium, which we calculate from the simulation in post-processing. Our non-equilibrium simulations are able to reproduce the relation between CO and H$$_2$$ column densities, and the relation between $$X_{\rm{{\small CO}}}$$ and metallicity, seen within observations of the Milky Way. We also compare to the xCOLD GASS survey, and find agreement with their data to our predicted CO luminosities at fixed star formation rate. We also find the multivariate function used by xCOLD GASS overpredicts the H$$_2$$ mass for our simulations, motivating us to suggest an alternative multivariate function of our fitting, though we caution that this fitting is uncertain due to the limited range of galaxy conditions covered by our simulations. We also find that the non-equilibrium chemistry has little effect on the conversion factor (<5 per cent) for our high-mass galaxies, though still affects the H$$_2$$ mass and $$L_{\rm{{\small CO}}}$$ by $$\approx$$25  per cent. 

Abstract The unprecedented infrared spectroscopic capabilities of JWST have provided high-quality interstellar medium metallicity measurements and enabled characterization of the gas-phase mass–metallicity relation (MZR) for galaxies atz≳ 5 for the first time. We analyze the gas-phase MZR and its evolution in a high-redshift suite of FIRE-2 cosmological zoom-in simulations atz= 5–12 and for stellar massesM_*∼ 10⁶–10¹⁰M_⊙. These simulations implement a multichannel stellar feedback model and produce broadly realistic galaxy properties, including when evolved toz= 0. The simulations predict very weak redshift evolution of the MZR over the redshift range studied, with the normalization of the MZR increasing by less than 0.01 dex as redshift decreases fromz= 12 toz= 5. The median MZR in the simulations is well approximated as a constant power-law relation across this redshift range given by $\log \left(Z/Z_{\odot }\right)=0.37\log \left(M_{*}/{\mathrm{M}}_{\odot }\right)-4.3$. We find good agreement between our best-fit model and recent observations made by JWST at high redshift. The weak evolution of the MZR atz> 5 contrasts with the evolution atz≲ 3, where increasing normalization of the MZR with decreasing redshift is observed and predicted by most models. The FIRE-2 simulations predict increasing scatter in the gas-phase MZR with decreasing stellar mass, in qualitative agreement with some observations. 

ABSTRACT Recent observations with JWST have uncovered unexpectedly high cosmic star formation activity in the early Universe, mere hundreds of millions of years after the big bang. These observations are often understood to reflect an evolutionary shift in star formation efficiency (SFE) caused by changing galactic conditions during these early epochs. We present FIREbox$$^{\it HR}$$, a high-resolution, cosmological hydrodynamical simulation from the Feedback in Realistic Environments (FIRE) project, which offers insights into the SFE of galaxies during the first billion years of cosmic time. FIREbox$$^{\it HR}$$ re-simulates the cosmic volume ($L=22.1$ cMpc) of the original FIREbox run with eight times higher mass resolution ($$m_{\rm b}\sim {}7800\, M_\odot$$), but with identical physics, down to $$z\sim {}6$$. FIREbox$$^{\it HR}$$ predicts ultraviolet (UV) luminosity functions in good agreement with available observational data. The simulation also successfully reproduces the observed cosmic UV luminosity density at $$z\sim {}6{\!-\!}14$$, demonstrating that relatively high star formation activity in the early Universe is a natural outcome of the baryonic processes encoded in the FIRE-2 model. According to FIREbox$$^{\it HR}$$, the SFE–halo mass relation for intermediate mass haloes ($$M_{\rm halo}\sim {}10^9{\!-\!}10^{11}\, {\rm M}_\odot$$) does not significantly evolve with redshift and is only weakly mass-dependent. These properties of the SFE–halo mass relation lead to a larger contribution from lower mass haloes at higher z, driving the gradual evolution of the observed cosmic UV luminosity density. A theoretical model based on the SFE–halo mass relation inferred from FIREbox$$^{\it HR}$$ allows us to explore implications for galaxy evolution. Future observations of UV faint galaxies at $$z\gt 12$$ will provide an opportunity to further test these predictions and deepen our understanding of star formation during Cosmic Dawn. 

ABSTRACT Observations show a tight correlation between the stellar mass of galaxies and their gas-phase metallicity (MZR). This relation evolves with redshift, with higher redshift galaxies being characterized by lower metallicities. Understanding the physical origin of the slope and redshift evolution of the MZR may provide important insight into the physical processes underpinning it: star formation, feedback, and cosmological inflows. While theoretical models ascribe the shape of the MZR to the lower efficiency of galactic outflows in more massive galaxies, what drives its evolution remains an open question. In this letter, we analyse how the MZR evolves over z = 0–3, combining results from the FIREbox cosmological volume simulation with analytical models. Contrary to a frequent assertion in the literature, we find that the evolution of the gas fraction does not contribute significantly to the redshift evolution of the MZR. Instead, we show that the latter is driven by the redshift dependence of the inflow metallicity, outflow metallicity, and mass loading factor, whose relative importance depends on stellar mass. These findings also suggest that the evolution of the MZR is not explained by galaxies moving along a fixed surface in the space spanned by stellar mass, gas-phase metallicity, and star formation rate. 

ABSTRACT Understanding what shapes the cold gas component of galaxies, which both provides the fuel for star formation and is strongly affected by the subsequent stellar feedback, is a crucial step towards a better understanding of galaxy evolution. Here, we analyse the H i properties of a sample of 46 Milky Way halo-mass galaxies, drawn from cosmological simulations (EMP-Pathfinder and Firebox). This set of simulations comprises galaxies evolved self-consistently across cosmic time with different baryonic sub-grid physics: three different star formation models [constant star formation efficiency (SFE) with different star formation eligibility criteria, and an environmentally dependent, turbulence-based SFE] and two different feedback prescriptions, where only one sub-sample includes early stellar feedback. We use these simulations to assess the impact of different baryonic physics on the H i content of galaxies. We find that the galaxy-wide H i properties agree with each other and with observations. However, differences appear for small-scale properties. The thin H i discs observed in the local universe are only reproduced with a turbulence-dependent SFE and/or early stellar feedback. Furthermore, we find that the morphology of H i discs is particularly sensitive to the different physics models: galaxies simulated with a turbulence-based SFE have discs that are smoother and more rotationally symmetric, compared to those simulated with a constant SFE; galaxies simulated with early stellar feedback have more regular discs than supernova-feedback-only galaxies. We find that the rotational asymmetry of the H i discs depends most strongly on the underlying physics model, making this a promising observable for understanding the physics responsible for shaping the interstellar medium of galaxies.